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Abstract 

We present a family of spherically symmetric multi-horizon spacetimes with a vacuum 
dark fluid, associated with a time-dependent and spatially inhomogeneous cosmological 
term. The vacuum dark fluid is defined in a model-independent way by the symmetry 
of its stress-energy tensor, i.e., its invariance under Lorentz boosts in a distinguished 
spatial direction (p r = —p for spherical symmetry), which makes the dark fluid essentially 
anisotropic and allows its density to evolve. The related cosmological models belong to the 
Lemaitre class of models with anisotropic fluids and describe a universe with several scales 
of vacuum energy related to phase transitions during its evolution. The typical behavior 
of solutions and the number of spacetime horizons are determined by the number of 
vacuum scales. We study in detail a model with three vacuum scales: GUT, QCD and 
that responsible for the present accelerated expansion. The model parameters are fixed 
by the observational data and by analyticity and causality conditions. We find that our 
Universe has three horizons. During the first inflation the Universe enters a T-region 
which makes the expansion irreversible. After the second phase transition at the QCD 
scale the Universe enters an R-region, where for a long time its geometry remains almost 
pseudo-Euclidean. After crossing the third horizon related to the present vacuum density, 
the Universe should enter the next T-region with inevitable expansion. 

PACS numbers: 04.70.Bw, 04.20.Dw, 04.20.Gz, 98.80.Hw 



1. Introduction 

Astronomical observations give a compelling evidence for the existence of a dark energy dom- 
inating our Universe at above 73 % of its density and responsible for its accelerated expan- 
sion due to negative pressure, p = wp, w < —1/3 [1-9] with the best fit w = — 1 [10-15] 
which corresponds to the Einstein cosmological term Xg^ related to the de Sitter vacuum 
Tfj, v = SirGp^g^ (A = 87rGp( vac )). 

The well-known cosmological constant problem has two aspects: (i) Quantum field theory 
predicts for p( vac ) the Planck scale p( vac ) = Ppi = 5 x 10 93 g/cm 3 . Confronting this with the 
observational value, p( vac ) — 3.6 x 10~ 30 g/cm 3 = 1.4 x 10 _123 ppi, creates the Fine- Tuning 
Problem, (ii) The inflationary paradigm needs a large value of p( va c) at the beginning of the 
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Universe evolution, typically of the GUT scale p( vac ) — Pgut — 5 x 10 77 g/cm 3 ; the observations 
indicate its much smaller value while the Einstein equations require p( vac ) = const. 

A typical solution which can be found in the literature, is to put P( vac ) = for some reason 
and to introduce a dark energy of non- vacuum origin which mimics the cosmological constant A 
when necessary. A lot of theories and models have been developed to describe a dynamical dark 
energy (for a comprehensive review see [9,10]). The alternative to the cosmological constant 
provided by quintessence assumes the existence of a hypothetical component of matter content 
with wq [17]. Q-models based on a scalar field, rolling down its self-interaction potential 

[18], were tested using different methods with WMAP-CMB data. This gave the constraint 
wq < —0.7, with the best fit Wq = — 1 [19] which evidently corresponds to the cosmological 
constant A. 

Quartessence models describe a transition from a dust-dominated stage to a late-time in- 
flationary stage. The Chaplygin gas model, with the postulated equation of state p = —A/p, 
gives a flat FRW model interpolating between p = and a negative pressure at late times [20]. 
It can be obtained in the model of a superfluid Chaplygin gas with the potential V(4>*4>) = 
M (0*0/p + p/0*0) which gives p = — 4M 2 /p [21,22]. In the holographic dark energy approach 
(see for a recent review [23]) with an interaction between dark matter with p = and dark 
energy [24], quartessence can be recovered as an isotropic perfect fluid [25], but the perfect fluid 
for holographic dark energy was found to be classically unstable [26] . The generalized Chaply- 
gin gas (GCG) model, p = —A/p a [27], was introduced to overcome difficulties with satisfying 
the CMB constraints [28]. The observational constraint on the parameter a, < a < 0.2, 
implies a little difference between the GCG and A = const [29]. 

Quintom cosmology describes the dark energy in the framework of the brane-world cos- 
mology by introducing two scalar fields, one being a quintessence and the other a phantom 
(for a review see [30]). The multiple- lamb da cosmology [31] describes the Universe as a kind 
of multiverse [32] filled with phantom energy; it may describe the evolution as a sequence 
of transitions between different inflationary stages. It is based on a perfect isotropic fluid 
with a time-dependent equation of state [33], introduced phenomenologically [31] and describ- 
ing phantom-non-phantom transitions [34]. One more approach to creating different effective 
scales of vacuum energy density rests on curvature-nonlinear multidimensional gravity with at 
least two extra factor spaces whose scale factors behave as scalar fields in four- dimensional 
space-time [35]. 

According to observational data, dark energy is well described by the inflationary equation 
of state, p = —p. The stress-energy tensor has the form 

T% = P(vac)^, P(vac) = COnst. (1) 

It represents a de Sitter vacuum which generates the de Sitter geometry responsible for accel- 
erated expansion independently of specific properties of particular models of p( vac ) • 

The Standard Models of cosmology and particle physics suggest a series of phase transitions 
that have occurred in the course of the expansion and cooling history of the Universe [36]. A 
connection between particle physics and cosmology predicts the first inflation related to a phase 
transition at the GUT (Grand Unification) scale. The first inflation solves the key problems of 
the standard Big Bang model ( [37] and references therein) and has been confirmed by CMB 
observations [36]. The Standard Model of particle physics predicts another phase transition at 
the QCD (Quantum Chromodynamics) scale of 100 to 200 MeV ( [36] and references therein) 
which can be related to a second inflationary stage [38]: A quasi-stable QCD vacuum state can 
lead to a short period of inflation (7 e-foldings) consequently diluting the net baryon to photon 
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ratio to its presently observed value. The second inflationary stage is considered in the model 
of thermal inflation [39], with a duration of about 10-foldings. Arguments for a second inflation 
at the QCD stage also exist in an effective model of a QCD phase transition which displays 
a high degree of supercooling at a critical temperature of the order of 100 MeV, so that the 
Universe increases exponentially during the quark-hadron transition [40]. 

The aim of this paper is to present and study a family of cosmological solutions to the 
Einstein equations describing a vacuum-dominated universe which several scales of vacuum 
energy related to phase transitions in the course of its evolution. 

The gauge non-invariance of quantum cosmology leads to a connection between the gauge 
and the quantum spectrum of a certain physical quantity which can be specified in the frame- 
work of the minisuperspace model. There exists such a particular gauge in which the cosmo- 
logical constant A is quantized. Transitions between quantum levels of the operator A can 
be related to several stages in the Universe evolution with different values of vacuum density 

P(vac) [41]. 

At the classical level, the key point is the algebraic structure of the source term in the 
Einstein equations. In a model-independent approach, a vacuum is defined by the symmetry 
of its stress-energy tensor [42-44], as suggested by the Petrov classification for stress-energy 
tensors. The Einstein cosmological term corresponds to the maximally symmetric de Sitter 
vacuum with p( vac ) = const. 1 The Petrov classification of stress-energy tensors implies the 
existence of vacua whose symmetry is reduced as compared with (1), which allows the vacuum 
energy to become time-dependent and spatially inhomogeneous [44,47,48]. A relevant class of 
stress-energy tensors describes a vacuum dark fluid [49] specified by the inflationary equation 
of state p a = —p in only one or two distinguished spatial directions, so that the vacuum dark 
fluid is intrinsically anisotropic. 

In the spherically symmetric case, a cosmological vacuum is specified by T\ = T r r (p r = —p) 
[43,44]. The radial direction is distinguished by the cosmological expansion. Regular solutions 
with source terms specified by T\ = T r r , necessarily have a de Sitter center [47]. In the case of 
two vacuum scales, at the center and at infinity, the source terms evolve from Ag^ to \g^ v with 
A C A [44,47]. The cosmological models belong to the Lemaitre class of models with anisotropic 
pressures. They are asymptotically de Sitter at the early-time and late-time stages [50,51]. 

In this paper we study, in a general setting, spherically symmetric spacetimes with several 
vacuum scales. The relevant Lemaitre models involve several de Sitter (inflationary) stages 
in the Universe evolution. We study in detail a cosmological model with three basic vacuum 
scales: the GUT and QCD scales, and that responsible for the presently observed accelerated 
expansion, p\ = (8irG)~ 1 \. We introduce a phenomenological density profile describing vacuum 
decay at each stage by an exponential function, as is typical for a decay, and fix the decay rate 
by the conditions of analyticity and causality. This approach allows us to reveal certain general 
features of our Universe including the number of its spacetime horizons. 

The paper is organized as follows: In Sec. 2 we introduce spherically symmetric vacuum 

1 Quantum field theory in curved spacetime does not contain a unique specification for the quantum state 
of a system, and the symmetry of the vacuum expectation value of a stress-energy tensor does not always 
coincide with the symmetry of the background spacetime [45]. In the case of de Sitter space, the renormalized 
expectation value of ( T M „ ) for a scalar field with an arbitrary mass m and curvature coupling £ is proved to have 
a fixed point attractor behavior at late times ( [45] and references therein), approaching, depending on m and 
£, either the Bunch-Davies de Sitter- invariant vacuum or, in the massless minimally coupled case (m = £ = 0), 
the de Sitter-invariant Allcn-Folacci vacuum. The latter case is peculiar since the de Sitter-invariant two-point 
function is infrared-divergent, and the vacuum states free of this divergence are 0(4)-invariant Fock vacua; the 
vacuum energy density in the 0(4)-invariant case is not the same (larger) as in the de Sitter-invariant case [46]. 
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spacetimes. In Sec. 3 we show how the number of vacuum cales determines the general features 
of a spacetime, including the number of horizons. In Sec. 4 we describe a transition from the 
static reference frame to geodesic reference frames representing Lemaitre cosmologies. Sec. 5 
presents a Lemaitre cosmological model with a vacuum dark fluid, T\ = T r r . Sec. 6 describes 
a model with three vacuum scales, GUT, QCD and present-day vacuum density, with the 
parameters fixed by the observational data. In Sec. 7 we summarize and discuss the results. 



2. Vacuum energy in general spherically symmetric spacetimes 

A general model-independent approach based on the Petrov classification of stress-energy ten- 
sors (SETs) defines a vacuum by the symmetry properties of its SET [42-44]. The Einstein 
cosmological term corresponds to the de Sitter vacuum with the SET 

t; = P s;, P = - P . (2) 

The medium specified by (2) is interpreted as a vacuum due to the algebraic structure of 
its SET (2). It has an infinite set of comoving reference frames, so that an observer cannot 
in principle measure his/her velocity with respect to it [42], which is an intrinsic property of 
a vacuum [52]. The Einstein equations imply V^T^ = 0, which leads to p = const for the 
de Sitter vacuum (2). The maximum symmetry of the vacuum SET (2) can be reduced while 
keeping its vacuum identity [43], and this inevitably leads (due to V V T" = 0) to a dynamical 
vacuum energy [44]. The vacuum SET with a reduced symmetry, such that only one or two of 
its spatial eigenvalues coincide with the temporal eigenvalue, represents a vacuum dark fluid 
with the equation of state p a = —p in the distinguished direction(s) [49]. 

The general time-dependent spherically symmetric spacetime is described by the metric 

ds 2 = e^dt 2 - e 2a dR 2 - r 2 dtt 2 , dtt 2 = d0 2 + sin 2 8d<p 2 , (3) 

where a, 7 and r are functions of R and t. The Einstein equations with source terms whose 
algebraic structure is specified by [43] 

7? = t; ( Pr = -p) (4) 

admit a class of regular solutions with a de Sitter center, = A5^ as r = 0, where A = 8nGpo 
corresponds to a certain fundamental scale of symmetry breaking p = p( vac ) at r = [47, 48] . 

In a comoving reference frame, the vacuum dark fluid specified by (4) is presented by the 
SET 

T ^ = diag(p, p, -p ± , -p ± ), (5) 

where p r and p± are the radial and transversal pressures, respectively. 

It can be easily shown that under these general conditions we necessarily have p = p(r) 
and hence also pj_ = p±(r). To begin with, the tensor (5) is invariant under any coordinate 
transformations in the (t, R) 2D subspace, which is just a definitive property of a vacuum. 
Moreover, if there is no material source of gravity other than (5), the system satisfies all 
conditions of the generalized Birkhoff theorem [53,54], whence it follows that there exists a 
coordinate frame (t, R) in which the metric (3) is t-independent, and consequently p and p± 
are functions of R alone. Let us show, however, that it is unnecessary to assume that (5) is 
the only source of gravity: it is sufficient to require that it does not interact with other kinds 
of matter, and thus the conservation law V U T^ = holds. 
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Indeed, in this case we have (dots and primes stand for d/dt and d/dR, respectively) 



r 



p + 2-(p + p ± ) = 0, (6) 



r 



p' + 2-(p + P± ) = 0. (7) 

If r = r(R), hence r = 0, Eq. (6) immediately gives p = p(R). If, on the contrary, r is t- 
dependent, so that in the most general case we can suppose p = p(R,r), then (6) and (7) 
combined lead to dp/ OR = 0, as was asserted, and from (6) it then also follows 

f \ rd P (*\ 

P±=P±{r) = -p-- — . (8) 

From (4) it follows G$ = G\ due to the Einstein equations. In the static reference frame, 
using the Schwarzschild coordinate R — r, this gives a' + 7' = in (3), whence, choosing the 
appropriate time scaling, we get a + 7 = 0, and the metric (3) takes the form 

ds 2 = A{r)dt 2 --JTK- A (r) = e 2l{r) (9) 

where the metric function A(r) is given by 

2M(r) , s 

Ar = 1 — , 10 

r 

and M(r) is the mass function 

rr 

M(r) = 4vr / p(x)x 2 dx. (11) 
Jo 

We adopt the weak energy condition, i.e., a non- negative density for any observer on a time- 
like curve, which is natural for cosmological models describing the evolution of our Universe. 
This condition requires p± + p > which leads, by (8), to a monotonically decreasing density 
profile p(r) [47]. This fact, together with the number of vacuum scales at which p± satisfies 
p± = —p, determines the generic behavior of the metric function A(r), and as a result the 
maximum number of spacetime horizons. The actual number of horizons is determined by the 
specific form of the profile p(r). In Section 6 we will introduce a density profile appropriate 
for three vacuum scales and show how the number of horizons in our Universe follows from the 
observational constraints. 

Eqs. (4) and (8) give an r-dependent equation of state. It is evident that an anisotropic fluid 
needs two different equation-of-state parameters w a = p a / p. In our case, w r = p r /p = —1 due 
to Eq. (4), and w± = p±/p = — 1 — (r/2)d(\np)/dr due to Eq. (8). The parameter w± satisfies 
w± > — 1 since p(r) is a monotonically decreasing function. The parameter w± approaches 
w± = —1 as r 0, as r —7- 00 where p approaches the present-day vacuum density p\ = const, 
and also at each intermediate inflationary stage with p( vac ) = const. 

The mass function (11) can be related to the Schwarzschild mass if we separate in (11) 
the presently observed vacuum density p\ as a background density. It is possible because 
p(r) is monotonically decreasing function, and p\ is its minimum value. If we introduce 
p{r) = pd + Pa, where pd is a dynamical density decreasing smoothly from the value at 
the center p = (87rG)~ 1 A to zero at infinity, then the mass function (11) takes the form 
M(r) = An Jq pd(x)x 2 dx + r 3 A/6 and contains, in the limit r — > 00, the Schwarzschild mass 
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M = An J °° pd(r)r 2 dr (in the Schwarzschild geometry it is measured by the Kepler law in the 
Newton limit of the Schwarzschild solution at its asymptotically flat infinity). The mass func- 
tion (11) differs from the proper mass in curved spacetime, which is obtained by integration 
with the proper volume element dV = y^gc/ 3 £, where 3 g refers to the determinant of the 
spatial metric. The difference represents the binding energy [55], called also the gravitational 
mass defect [52]. The dynamics of Lemaitre class models is determined by the mass function 
(11), therefore we do not consider here the proper mass. Let us note, however, that in the case 
of a cosmology which is asymptotically de Sitter at infinity, the notion of the total conserved 
(proper) mass of the Universe can be introduced, because any asymptotically de Sitter space- 
time must have an asymptotic isometry generated by the Killing vector d/dt, and there exists 
the notion of a conserved total mass of the spacetime as computed at the future infinity [56]. 

3. The number of horizons and the number of scales 

A typical behavior of the metric function A(r) in (9) is dictated by the dynamics of the transver- 
sal pressure p± in its source (5), which determines the maximum number and nature of its 
extrema, determining in turn the maximum number of horizons. 

According to the Einstein equations, the transversal pressure p± can be expressed in terms 
of the metric function A(r) as 

8„Gp ± = lA"+-. (12) 
2 r 

At an extremum of A(r), A' = 0, hence, if p± > 0, this extremum is a minimum, and this 
minimum of A is unique in the domain where p± > (otherwise there would be a maximum 
between two minima). Assuming that p± is normally positive (so that the strong energy con- 
dition holds) and becomes negative only at distinguished spatial domains related to particular 
stages with a de Sitter vacuum behavior, we fix the number of zeros of p± and restrict the 
maximum number of zeros of A(r), i.e., of spacetime horizons. One vacuum scale is related to 
the de Sitter center where p± is negative and A(r) has a maximum. If there is no other such 
scale, the transversal pressure changes its sign once, and A(r) has one minimum at which A(r) 
can be negative. In the asymptotically flat case A(r) — > 1 as r — > oo, hence one zero of p± can 
result in two zeros of A(r), and the spacetime can have, as a maximum, two horizons [47]. If 
there is also a de Sitter asymptotic at infinity, this gives another domain where p_i_ is negative. 
Then A(r) has two maxima and can have one minimum in between where p± is positive. In this 
case p± changes its sign twice, and the single minimum of A(r) leads to at most 3 horizons [50]. 

Each intermediate vacuum scale produces two more zeros of p± and at most two horizons. 
Hence, for n vacuum scales with negative pressure we find at most 2n horizons in asymptotically 
flat spaces and at most In — 1 horizons in asymptotically de Sitter spaces. Note that it is the 
maximum number of horizons, and it will be smaller if, under the same behavior of p±, the 
metric function is A(r) > at some of its minima or A(r) < at some of its maxima. 

The dynamics of the transversal pressure also determines a typical behavior of the equation- 
of-state parameter w±. For example, in the case of three vacuum scales of interest here, Pgut 
at the center, Pqcd and p A , appearing successively, there can exist two domains where p±_ is 
positive. The parameter w± = p±/ p takes the value w±_ = — 1 during each inflationary stage, 
and it can be positive during the transitions Pgut — > Pqcd, and Pqcd _ > Pa- 

The requirements of regularity at the center and the dominant energy condition do not 
restrict the total number of horizons. This can be seen from the following example. 
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Let p > and consider two extreme equations of state compatible with the dominant energy 
condition: (a) p± = —p and (b) p± = p. In the case (a) the most general static, spherically 
symmetric solution is Schwarzschild-de Sitter; in the case (b) the SET structure is the same 
as for a radial electromagnetic field, so that we arrive at the Reissner-Nordstrom metric. The 
metric function A in these two cases is 

(a) A{r) = 1 - — - H 2 r 2 , (13) 

(b) A(r) = 1 - ^ + t, (14) 

with constant parameters m, p, H, q. (Note that here q is not a charge, the notation is adopted 
for convenience). 

Now, suppose that in the space-time with a regular de Sitter center the equation of state 
(a), leading to the exact de Sitter metric, holds in a finite interval of r until A reaches zero (that 
is, for r < hi = 1/H). Beyond this horizon let us take the equation of state (b), so that the 
metric function A has the form (14). It matches to the solution in the previous interval of r if 
the constants p and q are found from the continuity conditions for A and A' . Since A'(hi) < 0, 
this necessarily means that h\ is the inner horizon of the metric (14). With growing r it will 
eventually reach the outer horizon r = /i2 with A' > 0. At this point we again switch the 
equation of state to (a), so that the next interval will be described by A{r) given in Eq. (13), 
with m and H determined from the continuity of A and A'. At the inevitable next horizon 
we will again join the Reissner-Nordstrom metric in the same way and so on. 

The process can be continued as long as one wishes. Its feasibility is guaranteed by the 
following facts verified by a direct inspection: 

(i) Given any h > and C < 0, one can always find such p and q that the function (14) 
satisfies the conditions A(h) = and A'(h) = C. Thus a next Reissner-Nordstrom segment can 
always be joined at the points hi, h 3 , etc. 

(ii) Given any h > and C > such that hC < 1, one can always find such m and H that 
the function (13) satisfies the conditions A(h) = and A'(h) = C; moreover, the condition 
hC < 1 always holds if h is the greater of two zeros of the function (14) and C = A'(h) for the 
same function. Thus a next Schwarzschild-de Sitter segment can always be joined at the points 
h%, h±, etc. 

The process can be stopped at any stage. If the last equation of state is (a), we obtain an 
asymptotically de Sitter model with an odd number of horizons; on the contrary, (b) leads to 
an asymptotically flat model with an even number of horizons. 

The density profile is continuous but contains fractures (jumps of the derivative p'). Where 
p± = —p, we have p = const, while where p± = p, the function p(r) behaves as 1/r 4 . The 
fractures can, however, be smoothed by arbitrarily small additions to p(r) without changing the 
whole qualitative picture, which will then correspond to an entirely smooth density distribution. 

Each plateau in the density profile must, from a physical viewpoint, manifest an intermediate 
energy scale, eventually connected with some phase transition. We can anticipate that the 
existence of such scales can appreciably complicate the set of possible spacetime structures. 



4. Transition to Lemaitre reference frames 



Consider the general static metric (9) that solves the Einstein equations with the SET (5). A 
transition to the geodesic coordinates (R,t), where r is the proper time along a geodesic and 
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the radial coordinate R is the congruence parameter, different for different geodesies, can be 
described in a general form. A radial timelike geodesic in the metric (9) satisfies the equations 

fdr\ 2 , . . dt E 

(_) = E>-A(r), - = W) , (15) 

where the constant E is connected with the initial velocity of a particle moving along this 
particular geodesic at a given value of the congruence parameter R. In general, E = E(R), i.e., 
it is different for different geodesies. 

Eqs. (15) give two of the four components of the transition matrix \\d(t,r)/d(T, R)\\, namely, 
f and t (dots and primes stand for d/dr and d/dR, respectively) since this partial differentiation 
occurs along the geodesies: 



r = ±^E 2 (R) - A(r), t = E(R)/A{r). (16) 

A relation between the other two components, t' and r', can be found from the condition g T R = 
when we substitute dt = tdr + t'dR and dr = rdr + r'dR into the metric (9): 

_ y/E*(R)-A(r) 

f ~ E(R)A(r) 1 j 

It remains to determine r'(R, r), which can be done by using the integrability condition (d T d^ — 
dnd T )r = 0. The latter takes the form of a linear first-order differential equation with respect 
to r' = y(R, r): 

yd r A EE' 

dry = ~2iE 2 —A ) + W^A- (18) 



Solving it, we obtain 

dr 



y 



r'(R,r) = VE 2 -A 



fo(R) + EE' 



(E 2 - A{r)f /2 



(19) 



The other integrability condition (d T dn — dnd T )t = holds automatically if (18) holds. The 
functions t(R,r) and r(R,r) can now be found by further integration of Eqs. (16)— (19). The 
resulting metric can be written as follows: 

d s 2 = dr 2 - r ' {I }; T } 2 dR 2 - r 2 (R, r)d{l 2 . (20) 
E 2 (R) 

One can see how this procedure works using de Sitter space as an example. Its static form 
is (9) with A{r) = 1 — H 2 r 2 , H = const. We will choose three different families of geodesies 
such that 

E(R) = Vl- KR 2 , K = 0,±1, (21) 

and show that their corresponding reference frames represent the three well-known forms of the 
de Sitter metric as isotropic cosmologies with different signs of spatial curvature (see, e.g., [57]). 
Indeed, integrating the first relation in (16) as an equation for r = r(r, R) and properly choosing 
the arbitrary function of R that appears as an integration constant, we obtain the following 
expressions for r: 

r(R,r) = (R/H) x {cosh(ifr), e Hr , sinh(tfr)}, (22) 
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where the expressions in the curly brackets are ordered according to K = 1, 0, —1. Substituting 
them into (20), we obtain the metric in the form 

ds 2 = dr 2 - a\r) + R 2 <Kl 2 ) , 

a(r) = (l/H) x jcosh(iJr), e Hr , sinh(tfr)}, (23) 

as was intended. One can also verify that the expression (19) with E(R) given by (21) (provided 
the function fo(R) is chosen properly) coincides with the expression for r' obtained directly from 
(22) in all three variants. So the transition has been completed. 



5. Lemaitre cosmology with vacuum dark fluid T\ = T r r 

We have shown above that the behavior of the static metric function Air) is dictated by the 
number of vacuum scales, and that the static spherically symmetric metric (9) can always 
be transformed to the Lemaitre form. Hence the cosmological evolution in this case can be 
described by a model from the Lemaitre class satisfying the condition p r = —p that specifies a 
vacuum dark fluid, with the appropriate choice of the density profile p{r) modelling smoothed 
jumps between different values of p(vac)> which will be discussed in the next section. 
A Lemaitre class model is described by the line element [52] 

ds ? = dT 2 _ e 2v(R,r) dR 2 _ r 2( R) r)dQ2) (24) 

The coordinates R, r are the Lagrange (comoving) coordinates. The coordinate r measures 
the proper time along the world lines of a fluid. The function r(R, r) corresponds to the Euler 
coordinate which is called luminosity distance. 

For the metric (24), the Einstein equations with the SET (5) read [52] 

8nGp r = -J (e" 2 V 2 - 2rr - f 2 - l) , (25) 

SnG P± = — (r"-rV)-- — £> — i> a — (26) 

e~ 2y 1 
87rGp = — (2rr" + r' 2 - 2rrV) + — (2rrz> + r 2 + l) , (27) 

2e~ 2v 

SttGTJ' = Or' - r'v) = 0, (28) 

r 

where dots and primes stand for d/dr and d/dR. 

The component T[ of the SET vanishes in the comoving reference frame since there is no 
momentum in the radial direction, and Eq. (28) is integrated giving [52,58] 

^ ^ TTwr (29) 

where f(R) is an arbitrary function. Putting (29) into (25), we obtain the equation of motion 
r 2 + 2rr + 8irGp r r 2 = f(R). (30) 
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Taking into account that p r + p = 0, the first integration of (30) gives 

r = ™m +m+ m (31) 

where the mass function M(r) is defined by 

M(r) = 4vr f p(x)x 2 dx. (32) 



o 



The arbitrary function F(R) (an "integration constant" parametrized by R) should be 
chosen equal to zero for models regular at r = since M(r) — > as r -> where p(r) — >■ p < oo. 
The second integration of Eq. (30) gives 

t-t (R)= [ -. dT . (33) 

J y/2GM(r)/r + f(R) 



The new arbitrary function t (R) due to this integration is called the bang-time function 
For example, in the case of the Tolman-Bondi model for dust (p r — pj_ — 0), the evolution is 
described by r(R, r) = [9GM(R)/2] 1 ^ 3 [t — r (i?)] 2 / 3 , where t (R) is an arbitrary function of R 
representing the Big Bang singularity surface at which r(R,r) = [60]. 

In the presently considered regular case, asymptotically de Sitter in the R-region near r = 0, 
the evolution starts from the timelike regular surface r(R,r) = rj. For f(R) > 0, the bang 
surface is r(R, r) = 0, and the solution (33) near this surface reduces to 

t — To(R) = I — F =======, (34) 

J ^/r 2 /r£ + f(R) 

where 



r °=\te (35) 

is the curvature radius at r = 0, and po is the density at r = 0. In the case f(R) < 0, the 
bang surface is r(R,r) = where r& satisfies 2GM(r)/r + f{R) = 0. For small values of 
f(R) we can apply (34) which gives r = r ^y—f(R) cosh [(r — r (i?))/r ]. For f(R) > we get 
r = r oy /f(R) sinh[(r - T (R))/r ], and r = r exp[(r - r (i2))/r ] for = 0. 

Different points of the regular timelike bang surfaces start at different moments of the 
synchronous time r, so that the bangs are non-singular and non-simultaneous. 

For f(R) = (parabolic motion), Eq. (34) gives at small r the expansion law 

r = r e^ T ~ T ° m/ro (36) 

and 



e 2 » = - 



dr Q (Ry' 



dR 

The metric takes the FRW form with the de Sitter scale factor 



(37) 
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where the variable q = e T °( R " ro is introduced to transform the metric to the FRW form. In 
accordance with (36), it describes a non-singular non- simultaneous de Sitter bang from the 
surface r(r — t (R) — > — oo) = [61]. 

The inflationary stage is followed by an anisotropic Kasner-like stage: One scale factor, 
corresponding to the transversal direction, is given by r(R, r), and the other, corresponding to 
the radial direction, is proportional to r' according to (29), its particular form depending on 
the density profile p(r) and the choice of arbitrary functions of R. If we choose f(R) = and 
t(R) = R, the metric can be approximated by [50,61] 

ds 2 = dr 2 - (r + R)~ 2/3 K(R)dR 2 - L(r + i?) 4/ W, (39) 

where K(R) is a smooth regular function and L is a constant. 

A similar behavior can be found for a density profile which approximates phase transitions 
with several scales of vacuum energy (see below). At each transition, an inflationary stage is 
followed by an anisotropic Kasner-like stage. 

The generic behavior of the considered Lemaitre class solutions is related to the function 
A(r) expressed by (10) in terms of the mass function (32). Given the density profile p(r), its 
specific form entirely determines the detailed properties of A(r) and, in particular, the number 
of horizons. 



6. Lemaitre cosmology with GUT and QCD phase transitions 
6.1. Basic features 

According to the conventional scenario, the first inflationary stage corresponding to the GUT 
phase transition occurred at the GUT scale -Equt ~ 10 15 GeV, the relevant density being 
p GUT ~ 2.3 x 10 77 g/cm 3 . It was followed by a decay of vacuum energy resulting ultimately in 
a radiation-dominated stage. The next phase transition which could drive the second inflation 
[38-40] which occurred at the QCD scale -Eqcd ~ (100^-200) MeV, at about 10~ 5 seconds after 
the Big Bang, when the Hubble radius, du = c/H, was about 10 km [36]. The density pqcd is 
smaller by a factor of (£ , qcd/-E'gut) 4 than the GUT density pgut = Po, i.e., of the order of the 
nuclear matter density. The last inflationary stage corresponds to the presently observed dark 
energy density p\ which is about 107 orders of magnitude smaller than the GUT density. 
This situation can be modelled by the density profile 

P = Po [l - (1 - #i) exp(-r"/r") - (B 1 - B 3 ) exp(-r 3 7r")] , (40) 

where n, Bi, B 3 , n, r 3 are constants, for which we adopt: 

Bi = Pqcd/Po ~ 10" 64 , B 3 = Px / Po w 10~ 107 ; r < r x < r 3 . (41) 

The exponential function in (40) is chosen as a typical one for decay processes. The parameter 
n characterizing the decay rate will be fixed below by the conditions of analyticity and causality. 

We choose f(R) = in Eqs. (29)-(33) because in this case each 3-hypersurface r = const 
is flat, with zero curvature [62], which guarantees fulfilment of the spatial flatness condition 
Q = 1 required by the observational data. 

Under the above choice, the model undergoes the following stages: 
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(a) r <C n: the first inflation, p ps po; the mass function (32) is approximated by M{r) = -yr 3 , 

and Eq. (33) yields, in agreement with (36), 

r-To(R) ~r In-. (42) 

(b) r ~ ri: end of the first inflation since the second term in (40) becomes significant. 

(c) r!<r< r 3 , so that 

p^p (B 1 + r n 1 /r n ). (43) 
This stage in turn splits into two periods. As long as r is sufficiently small, 

r < r 2 , r 2 = nB^ 1/n = 10 64/ Vi, (44) 

the second term in (43) is dominant, so that p(r) rapidly decreases. It is an intermediate 
period between the first and second inflation. At r = r 2 the two terms coincide, and at 
r > r 2 we have p ~ Bip = Pqcd = const, which corresponds to the second inflation. 

(d) r ~ r 3 : end of the second inflation since the third term in (40) becomes significant. 

(e) r ^> r 3 : the density is 

p«p (5 3 + £ir 3 7r n ). (45) 

Similarly to stage (c), at some value of r, namely, at r = r 4 defined by 

r 4 = r^Bi/Bz) 1 /" = 10 43/n r 3 = 10 107/n r! (46) 

the two terms in (45) are equal. At r 3 < r < r 4 , we have one more intermediate period 
where p(r) rapidly decreases, while at r > r 4 it approaches a constant corresponding to 
the present-day dark energy density. 

The time elapsed between the first and the second inflation, r 2 —T\ (we denote Tj = r(rj)) is 
estimated by integrating between r\ and r 2 in Eq. (33). To this end, we find the mass function 
in the same interval: 



Mir) = Ani^J + j j p(r)r 2 dr. 



In the first term we take p ~ Po while in the second one, in accord with (43), we approximate 
p(r) by p (ri/r) n . Hence, 

M W ^ w 3 + l^ f i i y 



3 """i 1 „_ 3 V r »-3 r «-3 

In the last term, almost in the whole interval of interest, r 3> r 1; therefore for our estimation 
purpose we can neglect the last term thus obtaining a constant value of M, 

M(r) ^M 2 = ——— porl = const. (48) 
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Substituting it into (33), we obtain 



/4(n-3) /r 2 \ 3 / 2 ti . 

T2-t 1 *\ -±- -r - (49) 



9n Vri7 

The first restriction on the parameter n is evident: n > 3. The second constraint follows 
directly from the dominant energy condition which requires p± < p and guarantees that the 
speed of sound never exceeds the speed of light, thus maintaining causality in the course of 
evolution. A simple analysis of Eq. (8) for several vacuum scales shows that the difference p—p± 
as a function of r decreases at each transition starting from 2p with p = const. Let us introduce 
the function /dec = P ~ P± characterizing the dominant energy condition. According to (8), 
/dec = 2p + rp'/2. It should be a decreasing function since p(r) is monotonically decreasing, 
and its derivative p ! is negative. During the first transition, this function should decrease from 
2p = 2pcuT to Ipqcd- For the density profile (40) we have 

f(r) = Wo" 1 = 2 [l - (1 - SOe-f /»* (l + \ (£)")] (50) 

It should be non- negative and decreasing from 2 to 2Bi, where B\ = Pqcd/ Pgut is given by 
(41). For r ^> n the exponent in (50) can be presented as a series in (ri/r) n which gives 



f = 2B 1 + 2(1 - B,) 



M)GF)"+M)(7f + °(7 



^ 3n 



(51) 



The condition of non-negativity of this function is n < 4. The derivative is given by 



'1 - Bi) n fri\ n 



f'= ^ r^ -p) e - r " /r " n -4-n(-) (52) 



ri \ n 



The condition n < 4 guarantees a monotonic decrease of the function p — p± during the first 
transition. It is easy to show that this concerns also the second transition at which the function 
/dec(^) decreases from 2p QC D to 2p X - 

The condition n < 4 thus provides non-negativity of the function /dec = P ~P± during the 
cosmological evolution with the density profile (40). Therefore we fix n = 4 as the only integer 
compatible with analyticity and causality. 

Let us note that at the end of both transitions, for r ^> n and r ^> rs, the density in (40) 
behaves like p oc r -4 , in a way typical of radiation in FRW cosmology where pa 4 = const and 
also agrees with our qualitative analysis in the second part of Section 3. 

With n = 4 we get from (47) 

M(r) = \ Wl 3 + 4 W 4 (---)■ (53) 
3 \ri r/ 

The value of M 2 in (48) is now 

16 

M 2 = — Wi- (54) 



Substituting it into (49), we obtain 
3 W 



T2 - TlS5 r£f^ 3/2 . 10 - r „. (55) 
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Recalling that r ~ 10 s Z p i, we find r 2 — T\ ~ 10 32 / p i/c ~ 10~ n s. It is of interest that this 
estimate does not depend on the particular choice of the free parameter n or, equivalently, the 
number of e- foldings N e := ln(ri/ro). 

Furthermore, if the second inflation contains 7 e-foldings [38], it means that r 3 ~ 10 3 r 2 . It 
is then also easy to find r 4 , the value of r at which the DE density has reached its modern 
value, from (46): r 4 ~ 10 43//4 r 3 . To estimate the duration of the second inflation r 3 — r 2 and the 
time r 4 of the onset of the latest A-dominated stage, we should integrate in (33) from r 2 to r 3 
and then to r 4 . Acting in the same manner as in finding r 2 , we see that the main contribution 
to the mass function comes from the range r < r 2 and is given by (54), hence the duration of 
the second inflation is 



Lastly, for r 4 we obtain the same relation as (56) but with r 3 replaced by r 4 . It results in 

r 4 - r 3 ~ 10 10 s ~ 1000 years. (57) 

Let us note that all these times are practically independent of the number of e-foldings iV e 
during the first inflation. However, the duration of the later period up to the present epoch 
does depend on N e . Namely, since at r ^> r 4 we have p ~ -B 3 po = const, integration in (32) at 
large enough r gives 

M{r) « (4vr/3) J B 3 p r 3 , (58) 

and integration in (33) yields immediately 

r - t (R) « 10 18 s-ln(r/r,), (59) 

where r* « 10 9 r 4 « 10 36 ri is the value of r at which the contribution of r 3 to the mass function 
begins to exceed M 2 . 

The qualitative behavior of the vacuum density profile (40) is shown schematically in Fig. 1. 
p/p 

1 .0 




10 107 



O.O 

Figure 1: Typical behavior of the vacuum density during the Universe evolution. Here r\ is the end 
of the first inflation, r 2 and r 3 are the beginning and end of the second inflation, respectively, r 4 
corresponds to achieving the present day vacuum density p\. 

Eq. (59) gives the expansion law r = r if e <yT ~ T °^^ r * , and Eq. (29) gives for the second scale 
factor in (24) e 2u = (r 2 / 'rf)(dro(R) / dR) 2 . Introducing the variable q = e T °( R ^ r * , we transform 
the metric (24) to the FRW form 

rfs 2 = dT 2 _ r 2 e 2cr/r, ^2 + g 2^2) ? (6Q) 
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at the stage where the vacuum density achieves its present value. 

As we have seen, the only free parameter of the model is the number of e-foldings at first 
inflation N e = ln(ri/ro), where the characteristic de Sitter radius for the GUT scale vacuum is 
r ~ 2.4 x 10 -25 cm. The time interval corresponding to r\ is, according to (42), approximately 
10 -35 _/V e s. In the next subsection we evaluate an admissible interval for the parameter N e from 
the requirement of late-time homogeneity and isotropy. 

For the density profile (40) the transversal pressure is given by 



P± = Po 



l + (l-S 1 )(l + ^)e-^ + (S 1 -5 3 )(l + ^) 



2rt 



(61) 



It satisfies the equation p± = —p during each inflationary stage. It has two maxima, p± ~ 
0.213p at r m i ~ 1.2r 1 and p± ~ 0.213£>ip = 0.213pqcd at r m2 ~ 1.2r 2 and ultimately 
quickly achieves p± = —B 3 p = —p\. The behavior of the transversal pressure (61) is shown 
schematically in Fig. 2. 




Figure 2: Typical behavior of the transversal pressure during the Universe evolution. The quantities 
r\—r± are the same as in Fig. 1. 

The parameter w± = — 1 during the first inflation, then it rapidly increases to w± ~ 0.213 
at r m \ ~ 1.2ri, decreases to w± = —1, quickly increases to wj_ ~ 0.213 at r m2 — 1.2r 3 , and 
finally approaches w± = — 1 as p approaches p\. 

6.2. Late-time homogeneity and isotropy 

At r > r 4 the model evolution is governed by the small effective cosmological constant A = 
B3P0/ (SnG) and tends to a de Sitter regime, i.e., becomes homogeneous and isotropic. 

The degree of inhomogeneity can be characterized by the dimensionless parameter (r j p)dpj dr 
showing how the density p changes at a distance ~ r. By (45) with n = 4, at r 3> r 4 this 
parameter is approximately equal to 4r 4 /r 4 and rapidly decreases with growing r. 

Thus at r > r4 one can estimate the degree of anisotropy of our model as is conventionally 
done for homogeneous models, e.g., using the anisotropy parameter [63,64] 

A =^pib H ?- 1 ' ( 62 ) 

i=l 

where Hi = hi/ai are the directional Hubble parameters corresponding to the three scale factors 
ai(r), the dot stands for d/dr, and H = (Hi-\-H2+H^)/3 is the mean Hubble parameter (see [64] 
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for a discussion of different anisotropy characteristics). In our model with the metric (24), where 



„/2 



and f(R) = 0, these scale factors are a± 



is found from (33): 

r' = -y/2GM(r)/rT^(R). 
Using this, one obtains for the anisotropy parameter 
(M/M -3f/r) 2 r 3M 4 - Axp±r\ 



and a2 = a-3 = r. The expression for r' 



(63) 



A 



(M/M + 3r/r) 



3M 4 + 47rp 4 (2r 3 



(64) 



where M 4 = M(r 4 ) and p 4 = p(r 4 ) = 2B%po according to (45). At large r the parameter 



A 



but, as can be directly verified, this rapid decrease does not begin from r 4 but only 



from much larger values of r because at r ~ r 4 both the numerator and the denominator of 
(64) are dominated by the constant M 4 . 

To agree with CMB observations, the vacuum contribution must be already highly isotropic 
(A < 10~ 6 ) when r reaches the value of the scale factor r = r 5 ~ 10 25 cm corresponding to the 
recombination epoch with redshifts z ~ 1000. This requirement constrains the possible value 
of the free parameter N e = \n(ri/r ). Indeed, the condition A{r = r 5 ) < 10~ 6 gives 



3M 4 < 2V2nGp 



y 5 



10" 



(65) 



(taking into account that M 4 S> 7rp 4 r|). In turn, GM 4 is expressed in terms of r\. From (54) 
we know the value of Mo. = Mir?) = (16/3)7rpo^i- To find M 4 , we must integrate in (32) from 
T2 to r 4 ; it turns out, however, that this integration contributes only a relative correction of the 
order 10 -7 to Mo. Thus 



Ma 
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-np r 1 . 



Comparing (66) with (65), we obtain the constraint 

n < Sg /3 r 5 /10 w xl0" 37 r 5 = 3xl0" 12 cm, N e = ln(n/r ) < 30. 



(66) 



(67) 



Hence, after the recombination time corresponding to r = rs, the Lemaitre model (24) 
practically behaves as a homogeneous and isotropic FRW model (60). 



6.3. Evolution of the scale factors 

Now, having established the constraint (67), we can discuss the model evolution at all stages. 
For the spatially flat model, f(R) = 0, the line element (24) takes the form 

ds 2 = dr 2 - 6 2 (r, R)dR 2 - r 2 (r, R)dn 2 , (68) 

where we have introduced explicitly two scale factors: r(r, R) in accordance with (24) and 
b(r,R) = r'(r,R) in accordance with (29). For the integration "constant" in (34) we choose 
Tq(R) = —R to make the de Sitter asymptotics familiar. It is easily seen that in this case 
r( r + R) = dr/d(r + R)d(r + R)/dr = dr/d(r + R). 

Numerical integration of the Lemaitre equations during the first transition pgut - > Pqcd 
shows an exponential growth of both scale factors at the beginning when p± ~ p r = —p, 
followed by an anisotropic Kasner-like stage where the anisotropy of the pressures leads to an 
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r b 




Figure 3: The behavior of the scale factors during the first phase transition. The distances are 
normalized to 10 12 ro where r$ ~ 2.4 x 10 -25 cm is the characteristic GUT length scale for Mqut — 
10 15 GeV, and the time r + R is normalized to the GUT time tcuT = r /c ~ 0.8 x 10~ 35 s. 




Figure 4: The behavior of the velocities r(r + R), b{r + R) at the first transition. 

anisotropic expansion. The behaviors of the two scale factors r(r + R), b(r + R) during the first 
phase transition are shown in Figs 3 and 4. 

The behavior of the scale factors at the second phase transition is qualitatively quite the 
same for the density profile (40). 

The evolution of the scale factors during the whole Universe history is shown schematically 
in Fig. 5 plotted on the basis of Fig. 3 extended to the third inflationary stage (the hypersurface 
t + R = £4 in Fig. 5). The first inflation ends at the hypersurface t+R = t±, the second inflation 
occurs between r + R = £2 and r + R = £3. The sharp maximum in b(r + R) in Fig. 4, as well 
as the two maxima in the right panel of Fig. 5 are related to a maximum of p± seen in Fig. 2. 
At late times, due to isotropy, the evolution of the two scale factors is common and conforms 
to the standard flat de Sitter cosmology. 
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t, t 2 t 3 t 4 t=x + R 



Figure 5: Typical behavior of the scale factors r(r + R) and b(r + it!). 
6.4. Horizons 

We have described the Universe evolution from the viewpoint of the Lemaitre cosmological 
reference frame. Now we can find the number of horizons and present the global structure of 
spacetime. The proof in Sec. 3 concerned only the maximum number of horizons, and their 
smaller number is certainly possible, depending on the dispositions and durations of the phase 
transitions on the r scale. 

Let us find how the whole scenario looks from the static reference frame. The metric function 
A(r) calculated with the above fixed parameters is shown in Fig. 6, where the characteristic 
scales designated on the r axis are: 

T\ = 10~ 12 cm, N e = 29 at the end of the first inflation. 

r 2 = 1.2 x 10 4 cm at the beginning of the second inflation, 

r 3 = 1.2 x 10 7 cm at the end of the second inflation, 

r 4 = 6.7 x 10 17 cm at achieving the present-day vacuum (dark energy) density p\. 

The values of r$ = 10 25 cm and r§ = 10 28 cm approximately correspond to the recombination 
and to the beginning of the third (presently observed) inflation, respectively. The value r-j = 
7.5 x 10 29 cm corresponds to the cosmological horizon due to de Sitter vacuum with the density 
P\. The behavior of the metric function A(r), Fig. 6, testifies for the existence of three horizons 
for the case of the above fixed parameters corresponding to those of our Universe. 

The first horizon in our Universe, r_, is close to the de Sitter radius r$ corresponding to the 
GUT scale of the first phase transition. The second horizon r + ~ 7.4 x 10 13 cm distinguishes 
the additional essential length scale: for a long time our Universe evolves as the T + -region 
r_ < r < r + , i.e., the expansion is inevitable for all observers as dictated by the causal 
structure of spacetime. The present epoch gets into the the i?-region between two T + -regions. 
Near the point of achieving the present vacuum density the geometry generated by the vacuum 
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Figure 6: The metric function A(r) plotted with the parameters corresponding to our Universe 

fluid background becomes almost pseudo-Euclidean (A(r) ~ 1 — 10 -6 ). 

The global structure of spherically symmetric spacetime with three horizons, asymptotically 
de Sitter in the center and at infinity, is shown in Fig. 7 [50,65]. This picture shows how the 




Figure 7: Global structure of spherical space-time with three horizons 

manifold of events is seen by different observers. Let us note that the Carter- Penrose diagram in 
Fig. 7 covers the whole plane except for the squares bounded by the lines r = and r = oo. The 
lines r_ correspond to cosmological horizons for static observers (observers in hats) in the PL- 
regions < r < r_ denoted as TZC; r + are the black (white) hole horizons for static observers in 
the R-regions r + < r < r ++ denoted as U, and r ++ are their cosmological horizons. T + -regions 
r ++ < r < oo denoted as CC correspond to regular homogeneous anisotropic cosmological 
T-models of Kantowski-Sachs type [50,51]. The Lemaitre cosmological model shown in Fig. 7 
starts its evolution in the R-region 1ZC\ and goes consequently through the T + -region r_ < r < 
r+, the R-region r + < r < r ++ and the T + -region r ++ < r < oo. 



20 



7. Summary and discussion 

We have presented a general analysis for the case of several scales of vacuum energy p( vac ) 
corresponding to phase transitions involving inflationary stages. In our approach, the vacuum 
dark energy is described by a vacuum dark fluid defined by symmetry of its stress-energy tensor. 
In the spherically symmetric case, it is invariant under radial Lorentz boosts, acquiring the 
maximum symmetry of de Sitter vacuum only at inflationary stages. This makes the vacuum 
density P( vac ) time-dependent and spatially inhomogeneous. Cosmological solutions generated 
by the vacuum dark fluid belong to the Lemaitre class models with anisotropic pressures (their 
anisotropy follows directly from the variability of P( vac ))- 

The intrinsic properties of de Sitter space-time are responsible for an accelerated expansion, 
independently of particular properties of particular models of vacuum density associated with 
the cosmological constant. In a similar way, the intrinsic properties of geometries generated by 
a vacuum dark fluid can be in principle responsible for a variable vacuum density and make 
it possible to describe, on a common ground, the first inflationary expansion, the presently 
observed accelerated expansion, as well as inflationary stages related to phase transitions in the 
universe evolution, predicted by the Standard Model. 

To our knowledge, such an approach is applied for the first time for an analysis of the 
cosmological evolution in the case of several vacuum scales. 

The dynamics of cosmological models with a vacuum dark fluid is dictated by the number of 
vacuum scales: their number determines the behavior of the transversal pressure p±, which in 
turn determines the maximum number of horizons; their actual number depends on the model 
parameters. 

We have studied in detail the cosmological model for the case of three vacuum scales: 
GUT, QCD and that responsible for the presently observed accelerated expansion. We used 
a phenomenological density profile with a typical behavior for a cosmological scenario with 
inflationary stages followed by a decay of vacuum energy, which we describe by an exponential 
function typical of decay processes. The parameter characterizing the decay rate is tightly fixed 
by the requirements of analyticity and causality. Other parameters of the model are fixed by 
the values Pgut, Pqcd, Pa and = 1. The only free parameter is the number of e- foldings in the 
first inflation, which is estimated using the observational constraint on the CMB anisotropy. 

This model reveals the following features of our Universe: 

(i) Our spacetime has three horizons. 

(ii) The cosmological evolution starts with a non-simultaneous timelike de Sitter bang followed 
by a short stage of inhomogeneous and anisotropic expansion. 

(iii) During the first inflationary stage, the Universe quickly enters a T + -region which makes 
the expansion irreversible. The second phase transition occurs during this period. 

(iv) The Universe enters an i?-region (in which we actually live) when the second inflation had 
already terminated. Soon after that the vacuum density reaches its present value (it occurs 
at r rs r4 ~ 10 18 cm and r ~ T4 ~ 1000 years). For a long time the Universe geometry 
remains almost pseudo-Euclidean up to the scale factor r of approximately 3 x 10 27 cm, which 
corresponds to an age of about 3 x 10 9 years, when, according to the observational data, the 
present vacuum density p\ begins to dominate and the third inflation starts. Let us note that 
in a model taking into account matter and radiation, the vacuum density could achieve its 
present value later. 
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(v) After crossing the third horizon related to the present vacuum density (r^ ~ 7.5 x 10 29 cm), 
the Universe enters the second T + -region with inevitable expansion. 

We have seen that even our purely vacuum model can fairly well conform to the basic 
observational features of our Universe, which proves the viability of our approach. We hope 
that inclusion of matter and radiation in a further development of this approach will result in 
models more completely describing the observed cosmological picture. 

A general exact solution for homogeneous T-models has been found for a mixture of vacuum 
dark fluid with T\ = T T r and dust-like matter [51]. It presents the class of T-models specified by 
the density profile of a vacuum fluid. The solution contains one arbitrary integration constant 
related to the dust density. Numerical estimates for a particular model illustrate the ability of 
such models to satisfy the observational constraints [51]. A similar solution for Lemaitre class 
models is now under consideration. 
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